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Abstract. 1 present in this paper a method to subtract the bias due to source photon noise from visibilities 
measured with a single-mode optical interferometer. Properties of the processed noise are demonstrated and 
examples of subtraction on real data are presented. 
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1. Introduction 

The properties of source photon noise are well known. It 
follows a Poisson distribution whose variance is equal to 
the average total number of photons. In frequency space, it 
is a white noise with a flat average power spectral density. 
In most practical cases where the observables are linearly 
linked to the number of photons detected, photon noise 
can be directly averaged out from the data to reduce its 
variance. For some applications for which the observables 
are quadratically linked to the number of photons, the 
data suffer both from photon noise and from a bias linked 
to the variance of the noise. This is for example the case 
in speckle imaging techniques where the source spatial in- 
tensity distribution is recovered from the power spectral 
density of a short time exposure (Thiebaut (1994)). In as- 
tronomical optical interferometry, the observables are the 
modulus of the visibility and its phase usually expressed 
as a closure phase quantity. The visibility modulus can be 
obtained by integrating the modulus of the spectrum of 
interferograms. However, this estimator is biased by the 
power spectral density of noises as these add to the power 
spectral density of the fringe signal. An unbiased estima- 
tor of the modulus of the visibility is obtained by forming 
the squared modulus of the visibility as the power spectral 
densities of the noises can be independently estimated and 
subtracted. An example in interferometry is the computa- 
tion of the fringe squared visibility in the ABCD method 
where the white light fringe is sampled at four A/4 spaced 
optical path differences. An unbiased single fringe ABCD 
estimator is obtained by subtracting the source photon 
noise and the detector noise variances (when the noise 
has a flat spectrum the power spectral density is constant 



and equal to the variance) from the fringe power spectral 
density ( [Tango fc Twiss (1980)1 ). 

In single-mode interferometers, beams are spatially fil- 
tered by single-mode waveguides trading phase fluctua- 
tions against intensity fluctuations. A fraction of intensi- 
ties collected by each aperture can be measured to renor- 
malize interferograms to eliminate the fluctuations due to 
turbulence. The visibility estimator is no longer directly 
linked to the power spectral density of the fringe signal. 
I demonstrate in the following sections that the classical 
method (explaine d further in the paper) established by 
[Goodman (1985) can be rigourously extended to such ra- 
tios of physical noisy signals under certain assumptions to 
provide unbiased visibility estimators. Real data reduction 
cases are presented to illustrate the method. 

2. Principles of photon noise bias subtraction 

2.1. Description of the signals and assumptions 



I refer the reader to Coudc du Forcsto et al. (1997) for a 
full description of the principle to measure fringe ampli- 
tudes (also called coherence factors) with a single-mode 
fiber interferometer and for a detailed description of the 
fringe signal for coaxial interferometers. Vi sibility calibra - 
tion will be addressed in a separate paper ( Perrin (2002) ). 
I assume here, for sake of simplicity, that no calibration is 
required and the coherence factor is directly equal to the 
visibility. Here I will use the more general expression of 
the interferogram for a two-telescope interferometer: 



lix) = Pa{x) + Pb{x) + 2^PA{x)PB(x)m{x) 



(1) 
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where m is an oscillating function containing the fringes. 
Pa and Pb are the intensities coupled in the single- 
mode waveguide at each telescope and are called the 
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Fig. 1. Examples of raw signals acquired for Mira in October 2000 with FLUOR, a) and b) The two photometric 
signals PA{t) and Psit). c) Interferometric signal i{t). d) Power spectral density of the interferometric signal. Energy 
units are photon counts. The power spectral density is expressed in squared photon counts. The fringe peak is located 
between 350 and 400 Hz. Other peaks are due to the noise and disappear when averaging the power spectral densities. 
The source is resolved hence the fringe amplitude is small and hardly shows up in the interferogram signal. The 
fringe peak in the power spectral density is proportional to the squared modulus of the visibility which is obtained by 
integrating the peak. 



photometric signals, i is a function of the optical path 

difference or of time for coaxial beamcombincrs. For 
multiaxial interferometers, a; is a focal plane coordinate. 
In the following, I will deal with coaxial interferometers 
only and I will use time t as the variable. The method 
can be easily adapted to multiaxial interferometers. The 
photometric signals vary in time with turbulence and the 
modulation m varies with turbulence and with the optical 
difference which is varied linearly with time, coding the 
fringe signal in frequency space, i is measured in photon 
counts and is defined as the average signal one would 
obtain if no noise were present. 

In the photometric calibration method, the photo- 
metric signals need to be estimated. The estimated 
signals are filtered, the filtering function being adjusted 
to reject most of the noise and keep the intensity fluctua- 
tions due to turbulence only. Turbulent fluctuations are 
low frequency fluctuations (limited to frequency ranges 
of a few tens of Hertz) . I note Pa and Pb the estimated 
photometric signals suitably low-pass filtered. The impor- 
tant property of the filtered photometric signals is that 



they contain no energy at the fringe frequency and above. 
I call g the following gain function: 



9{t) = 



1 



2^PA{t)PB{t) 

I define the normalized interferogram: 



(2) 



in{t) = g{t).i{t) 

PA{t) + PB{t) 



(3) 



.JPA{t)PB{t) 



2^PA{t)PB{t) ^JPA{t)PB{t) 



m{t) 



The first term is mainly a low frequency signal whereas 
the second term is the high frequency signal containing the 
fringe modulation. I introduce the continuum function: 



PA{t) + PB{t) 



(4) 



2jPA{t)PB{t) 



The continuum function is the low frequency part of the 
normalized interferogram. It is the ratio of the arithmetic 
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Fig. 2. a) Gain function g{t). b) Continuum function c{t). c) Power spectral density of the gain function, d) Power 
spectral density of the continuum function. The gain function is homogeneous to the reciprocal of a number of photon 
counts. The continuum function has no unit. 



and geometric means of the photometric signals. If the two 
photometric signals are equal then the continuum function 
is equal to 1. It departs all the more from 1 as the pho- 
tometric beams get unbalanced. The visibility estimate is 
computed from the corrected interferogram: 



«„(i) - c{t) 

[PA{t) -Pl{t)] 



(5) 



[PBit)-PB{t) 



2jPA{t)PB{t) 



^fpimW) 



m{t) 



PA{t)PB{t) 



In the corrected interferogram the low frequency compo- 
nents due to turbulence are eliminated. If the photometric 
signals are perfectly estimated then the average value of 
the corrected interferogram is equal to zero and the os- 
cillating signal is properly renormalized. In the ideal case 
where data are noiseless, the corrected interferogram is 
equal to the oscillating function m which is proportional 
to the visibility. Although the filtered photometric signals 
do not contain energy at the fringes frequency, the gain 
and continuum function may contain some as the residual 
noise may have been redistributed in the frequency do- 
main by the non-linear combinations of the filtered pho- 
tometric signals. I make the assumption that this high fre- 



quency noise is negligeable by several orders of magnitude 
compared to the energy of the fringes. This assumption 
will be validated in Sect. 3 on real data. As a consequence 
of this assumption, z„ and icor share the same photon 
noise. This assumption is crucial for the success of the 
method. Finding an analytical solution to the problem of 
bias subtraction without this assumption is certainly a big 
issue. 



2.2. Method 

The method to subtract the photon noise bias is a direct 
generalization of that proposed by Goodman (1985) 
and I will use the same notations. In the following 
only photon noise is considered. Methods to subtract 
additive noises are well esta blished (see for example 
[Coude du Foresto et al. (1997) ). Detector and source 
photon noises being independent the detector noise 
could be added in the following derivation leading to the 
classical result on power spectral density bias by detector 
noise variance. The primary scope of this paper being 
source photon noise I have decided to keep equations as 
light as possible and not include detector noise in the 
equations. Detector noise will nevertheless be considered 
in the last section on real data. 



4 



G. Perrin: Subtracting the photon noise bias from single-mode optical interferometer visibilities 



0,8 




20 40 60 80 100 120 200 400 600 800 1000 

Time (ms) Frequency (Hz) 



Fig. 3. The corrected interferogram and its power spectral density. The corrected interferogram being the ratio of 
photon counts has no unit. 



I call i{t) a representation of the interferogram in 
which photon events are represented by Dirac functions. 

Assuming that the total number of photons in the 
interferogram is N for a given realization, then I can 
write: 



N 



(6) 



There arc two random variables in this expression: the 
individual photon detection times tk and the total number 
of photons N which is equal to N in average. A model of 
the normalized interferogram can be built: 



N 



N 



^n{t)=Y.9{t)S{t-tk) = J29itk)S{t-tk) 
fc=i fc=i 

Its Fourier transform is therefore: 



(7) 



N 



(8) 



fc=i 



Thus, the average spectrum of the normalized interfero- 
gram is therefore: 



N 



k=l 

N 



(9) 
(10) 



fc=i 



In the above expression, the average on tk does not depend 
upon k and the average normalized interferogram can be 
written: 



(/„(/)) = (TV (5(i)e-^"*0t)iv 



AT 



The statistics of the photon arrival time t are described 
by the probability density function ^ with i{t) the aver- 
age photon flux. Hence the average on the arrival times is 
equal to: 



1 



= ^/(/)*G(/) 



(12) 

(13) 



The -k symbol indicates a convolution. Capital letters are 
used for Fourier transforms. Replacing the average on t 
by this expression, I obtain for the average normalized 
interferogram: 



1 

G(/)*/(/) 



(/„(/)) = (iV)^-G(/)*J(/) 



(14) 
(15) 



As a clicck, the expression of the average interferogram 
after applying the inverse Fourier transform to the above 
expression is: 



{in{t)) = g{t).i{t) 



(16) 



Let us consider 7„ '{f) the power spectral density of the 
normalized interferogram with In\.f) = |-^n(/)P- In the 
ideal case of non noisy data, the integral of the power 
spectral density is proportional to the squared visibility. 
The integral of the average power spectral density is: 



N N 



J2J29{tk)9{ti)e- 



■2iTT{tk-tl)f 



(17) 



\fe=i 1=1 



An important assumption on photon events to derive the 
Poisson statistics is that they are not correlated. Arrival 
times tk and ti are therefore not correlated as long as 
k ^ I. This important property is used to split the above 
expression in two terms: 



(11) {iL'\f)) = 
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Fig. 4. a) Average corrected interferograms power spectral density, b) Average processed dark current scans power 
spectral density, c) Difference of the two to suppress the bias due to detector noise, d) Difference of the two with 
photon noise bias subtracted. The energy peak at low frequency in the processed dark signal power spectral density 
is due to the multiplication of the dark signals by the gain function. This peak causes a drop when the dark signal 
power spectral density is subtracted from that of the fringe signal. 



JV 



TV 



N 



J^iait.)? + 2Y.j:9itMi)e-'"^"'-''^n ) (18) 



(23) 



\fc=i 



fc=l Kk 



The first average is equal to: 



N 



N 



\fe=i 



tk' N 



\k=l 



(19) 



N 
,l{t) 



E/ mf-w'^n (20) 



[g{t)fi{t)dt 



(22) 



where I have used once again the probability density of 
the photon statistics ^ . The second average of Eq. (|l8|) 
can be written as the sum of factors yielding: 



N 



k=l Kk 



fe=l Kk 



N 



The averages on ti and tk can be substituted by the ex- 
pression of Eq. (nsh yielding: 



N 



\ \ k=l Kk 
I N 

2EEi^w)*G(/)i= 



tk.ti' JV 



k=l Kk 



(24) 



and: 



^ 1 

2EE]^w)*G(/)i' 



fe=l Kk 

N{N - 1) 



N 



N 



1 



|/(/)*G(/)p (25) 



-^|/(/)*G(/)P 

|/(/) ★G'(/)p is the power spectral density of the normal- 
ized interferogram. From the equation above I therefore 
derive an unbiased estimate of this quantity: 

\I{f)*G{f)\' = (nPif)) / [9{t)M)dt (26) 
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Since the normalized and the corrected interferograms 
share the same noise bias, the integral term of the above 
equation is also the bias of the power spectral density of 
the corrected interferogram. 

2.3. Comments 

The expression of the photon noise bias is very intuitive. 
i{t) is the average number of photons detected at time t. 
Being a Poisson statistics, it is also the variance of the pho- 
ton noise. When the signal is multiplied by the gain <?(<), 
the variance at time t becomes [g{t)]'^i{t). Photon events 
at different times being uncorrelated, the total variance of 
the photon noise is therefore equal to the integral of the 
local variance. 

Equation ( p6| ) demonstrates that the noise of the corrected 
interferogram remains a white noise whose mean power 
spectral density is constant. This property is also the re- 
sult of the independence of photon events. 
With the above result, the computation of the unbiased es- 
timate of the visibility is easy. The bias is simply obtained 
by co-adding the individual counts of the normalized in- 
terferogram. 

3. Example of bias subtraction on real signals 

The purpose of this section is to illustrate the theoretical 
results of the previous section and to show that the 
assumptions on the noise of the continuum and gain 
functions are correct. 

I have selected a series of scans of Mira observed 
with FLUOR in October 2000. The observations were 
carried out with a baseline long enough that the fringe 
contrast is of a few percent only. This is an interesting 
case because the source is very bright and the fringe 
contrast is small, hence the photon noise bias is relatively 
important and accounts for a large fraction of the visibil- 
ity if not corrected. 

The raw photometric and interferometric signals are 
displayed on Fig. ^. The intensities are in photon counts. 
The fringes are hardly visible and their amplitude is 
smaller than that of the photometric fluctuations. The 
fringes peak is visible in the power spectral density 
between 350 and 400 Hz. In Fig. ^ are presented the gain 
and continuum functions as well as their power spectral 
densities on a log scale. The noise level of the continuum 
at the fringes frequency range can be directly compared 
to that of the corrected interferogram in Fig. ^. A residual 
photometric fluctuation is visible on the corrected inter- 
ferogram hence the low frequency component. In the data 
reduction procedure, the detector dark current signals 
undergo the same correction process as the source signals. 
In case the dark currents would combine both detector 
dark signal and sky background if measured with a 
chopping technique, the dark signals would contain both 
background photon noise and detector noise. But both 



noises can be treated as additive noise, in particular the 
background photon noise does not need to be processed 
like the source photon noise. In the case of the FLUOR 
signals, background photon noise is totally negligible and 
the dark signals only contain the detector contribution. 
In Fig. ^ I have represented the average of power spectral 
densities of both the corrected detector dark current and 
corrected interferograms to reduce the noise on the power 
spectral densities. On the same figure, the difference of 
the two is plotted. At this stage, the bias due to source 
photon noise is not removed to make its magnitude 
obvious to the reader. The fringe squared visibility being 
obtained by integrating the power spectral density at 
the fringe frequency, the non-zero mean level under the 
fringe peak causes the bias. This is the bias due to source 
photon noise. This graph also shows an increasing residual 
noise at high frequency due to detector instabilities. The 
value of the bias has been computed with the method 
presented in this paper and subtracted as illustrated by 
the last graph of this figure. The power spectral density 
background level is now equal to zero in average showing 
that the bias has been correctly subtracted. 

4. Conclusion 

A method to subtract the photon noise bias from vis- 
ibility data acquired with a single-mode interferometer 
has been presented. An analytical expression for the bias 
can be established under verified assumptions. Other non- 
analytical methods based on the fit of the average level of 
power spectral densities may suffer from confusion with 
fringe signal and from detector instabilities and the ana- 
lytical method should be preferred. 
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